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As one test of the validity of the staggered-fermion fourth-root determinant trick, we examine the suppression 
of the topological susceptibility of the QCD vacuum in the limit of small quark mass. The suppression is sensitive 
to the number of light sea quark flavors. Our study is done in the presence of 2+1 flavors of dynamical quarks in 
the improved staggered fermion formulation. Variance-reduction techniques provide better control of statistical 
errors. New results from staggered chiral perturbation theory account for taste-breaking effects in the low-quark 
mass behavior of the susceptibility, thereby reducing scaling violations from this source. Measurements over a 
range of quark masses at two lattice spacings permit a rough continuum extrapolation to remove the remaining 
lattice artifacts. The results are consistent with chiral perturbation theory with the correct flavor counting. 
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1. INTRODUCTION 

The continuum Euclidean partition function 
for QCD, restricted to topological charge sector 
v with n/ degenerate flavors is 



Z v = J [dU] exp[5({7)] det[M(U)] n * 
det[M(U)} = m^Y[(\l+m 2 ) (1) 

n 

where U is the gauge field, S(U), the Yang-Mills 
action, and A„ are the nonzero eigenvalues of the 
Dirac operator D(U), where M{U) = D(U) + m 
pp. As the quark mass m vanishes, nonzero values 
of v are suppressed and the topological suscepti- 
bility 

X = (v 2 ) /V = Y, v2Z »l ZV 



for Eucildean four- volume V vanishes, 
perturbation theory predicts 



l/ X = 2/ M / 2 ^n. i /m l 



Chiral 



(2) 



♦Presented by C. DeTar. 



for several flavors of degeneracy rii and mass m,. 

The lattice partition function for staggered 
fcrmions has a similar form. But to compen- 
sate for the unwanted four-fold "taste" degen- 
eracy, it is common to replace det[M(L7)] n / by 
det[M(U)} n f /4 (fourth-root trick) with the hope 
that the continuum limit recovers Eq. Two 
principal lattice artifacts modify Eq. (J2J) 0- 
First, topological charge counting becomes am- 
biguous on the scale of the lattice spacing. Sec- 
ond, the spectrum of eigenvalues of the fermion 
matrix M(U) is modified so that would-be zero 
modes are shifted away from zero and the four- 
fold taste degeneracy required by the fourth-root 
trick is only approximate 
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The staggered fermion artifacts can now be ac- 
commodated within the framework of staggered 
chiral perturbation theory J3J, and, as we show 
from a numerical study with improved staggered 
fermions, improves the approach to the contin- 
uum limit, where we find a result consistent with 
Eq. ©. 

2. STAGGERED CHIRAL PERTURBA- 
TION THEORY 

We summarize the derivation of the leading 
S^PT formula for the susceptibility and refer the 
reader to Ref. 0] for details. We start from the 
Aubin and Bernard generalization of the formal- 
ism of Lee and Sharpe [S| to obtain the effec- 
tive chiral Lagrangian for staggered fermions to 
0(a 2 ): 

f 2 ii f 2 

C = L.Tr(d„U*d lt U)-!^1i[M(rf + U)] 



2 

j,2 
0/ 



The chiral field U describes N flavors of 4 tastes 
each and 4>qi is the flavor-singlet, taste-singlet 
field, M. is the taste-degenerate quark mass ma- 
trix, and V(U) is the dimension-six potential 
which breaks the continuum taste symmetry. 
Only one Goldstone pion in the sixteen-member 
taste-multiplet survives. The "pion" with the 
highest mass is a taste singlet. 



We work to quadratic order in the fields 
and in the mean-field approximation and follow 
Lcutwyler and Smilga yQ to obtain the suscepti- 
bility with a degeneracy of n, flavors of four tastes 
each: 



X 



/ 2 /16 



where the bare flavor-neutral, taste-singlet pseu- 
doscalar masses are 



m 



il 



2finii + a Aj 



in terms of the bare quark masses m, and the 
taste splitting Aj. To eliminate the taste degen- 
eracy, leaving only two degenerate u and d quarks 
and one s quark we put n u d = 1/2 and n s = 1/4, 
giving, finally 



X = 



f 2m l,i/ 



1 + mlr/2m 2 ss j + 3m 2 7 /2m 2 



(3) 



This leading-order formula is identical to that 
of continuum chiral perturbation theory, except 
that it requires the taste-singlet meson masses. 
It interpolates [Sj between the perturbative and 
quenched regimes and reproduces the Witten- 
Veneziano formula at infinite meson mass [Jj- 

3. SIMULATION AND ANALYSIS 

The topological charge density was measured 
as described in [2] on an ensemble of lattices 
(quenched and with up, down, and strange im- 
proved Asqtad staggered quarks) of two lattice 
spacings: "coarse" (a w 0.125 fm) and "fine" 
(a « 0.09 fm), as listed in the upper and lower 
tables respectively . Included in this analysis is 
a new fine-lattice 40 3 x 96 ensemble with u and d 
quark masses at 1/10 the nominal strange quark 
mass. 

To determine the susceptibility, we use a vari- 
ance reduction technique [Hj , based on the obser- 
vation that with large lattice volumes we gain sta- 
tistically by writing the susceptibility as a local, 
intensive observable, rather than a global observ- 
able: 
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Figure 1. Result of 0. The curves are from 
Eq.' © with (L.O. 2 + 1 + oo) and without (L.O. 
2+1) the mg term. 



where the measured topological charge density 
correlation function, is 

C mcas (r) = (p(O)p(r)) 

and p is the topological charge density operator. 
The correlator has the asymptotic form ^U] 

Cfit(r) = b n D(m n ,r) + b rr D(m n ,,r) + ... 

where D{m, r) is the Euclidean scalar propagator. 

We fit the correlator to this form for large |r| 
and use it to model the large distance behavior, 
thereby reducing the variance due to contribu- 
tions at large \r\: 
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4. RESULTS 



As we have seen, at lowest order in chiral per- 
turbation theory the staggered fermion artifact is 
removed by writing the susceptibility as a func- 
tion of the taste-singlet pion mass, rather than 
the Goldstone pion mass. The improvement can 
be seen by comparing Figs. and [21 

Residual lattice artifacts come from the lattice 
operator definition of topological charge density. 
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Figure 2. New result. 

It is plausible that they scale as 0(a 2 ). The con- 
tinuum extrapolation of Fig. [2]at fixed abscissa is 
based on that assumption. 
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